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Abstract 


This  investigation  considers  the  thermo -mechanical  effects  of  an 
asperity  traversing  at  a  high  speed  over  a  semi -infinite  medium  with  a 
thin  hard  coating  surface.  The  general  analytical  solutions  of  the 
mechanical  stress  state,  the  temperature  field  and  the  thermal  stress 
state  are  obtained  and  expressed  in  Fourier  transform  space.  The 
asperity  speed  is  sufficiently  high  such  that  the  thermal  stress  is 
much  larger  than  that  caused  by  the  mechanical  traction  of  asperity 
pressure  and  friction.  The  analysis  emphasizes  the  heating  effect  of 
the  friction  force,  which  leads  to  the  initiation  of  the  thermo¬ 
mechanical  cracking,  or  ^heat- checking"  in  the  surface  layer,  the 
substrate  or  their  interface.  For  hard  coating  layers,  the  initiation 

<  ’  *V»  r,  ■  ' 

of  a  crack  will  occur  at  a  depth  rf-y  .  (where  the  maximum  principal 

thermal  stresses  occur) .  The  depth  is  found  to  be  controlled 
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principally  by  the  Peclet  number.  It  is  found  that,  for  a  thick  ± 

coating  layer,  the  crack  will  occur  inside  the  surface  layer  where  the 
stress  state  is  little  affected  by  the  substrate,  for  which  the 
estimate  of  fracture  initiation  is  the  same  as  that  of  a  single 
material  of  the  coating.  Yet,  for  modified  surfaces  of  thin  coating 
layer  in  which  the  thickness  is  of  the  same  order  of  the  rj  ,  the 
interaction  between  the  surface  layer  and  the  substrate,  through  their 
differences  in  mechanical  and  thermal  properties,  greatly  influences 
the  combined  stress  state  near  the  interface,  and  thus  the  crack 
initiation  in  the  neighborhood  of  ^max-  The  investigation  employs  the 
Fourier  transform  technique  for  the  analytical  method  in  order  to 
facilitate  the  parametric  study.  The  complex  Fourier  transform 
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inversion  could  be  done,  within  the  strip  of  analyticity,  alor 
carefully  selected  path.  The  speed  of  numerical  integration  can 
optimized  with  the  proper  choice  of  the  offset,  c,  from  the  real 
which  is  the  theoretical  path  of  integration  for  Fourier  inverse 


tranform. 


Table  of  Contents 


Acknowledgements  . 

Abstract  . 

Table  of  Contents  . 

List  of  Figures  . 

Nomenclature  . 

Chapter  1  Introduction  and  Literature  Survey  . 

1.1  Problem  of  Thermo -mechanical  Cracking  in  Layered  Medium 

1.2  Parametric  Effects  . 

1.3  General  Theory  . 

Chapter  2  Problem  Statement  and  Solution  Technique  . 

2.1  Basic  Equations  . 

2 . 2  Governing  Equation  . 

2.3  Integral  Transform  Technique  and  General  Solution  . 

2.3.1  Fourier  Transform  Method  . 

2.3.2  Mechanical  Stress  Field  . 

2.3.3  Temperature  Field  . 

2.3.4  Thermal  Stress  Field  . 

Chapter  3  Numerical  Solution  . 

Chapter  4  Parametric  Study  . 

4.1  Asperity  Parameters  and  Critical  Depth  . 

4.2  Material  Parameters  . 

4.3  Coating  Layer  Thickness  Effects  . 

4.3.1  Effects  of  Differences  in  Mechanical  Impedance  ... 

4.3.2  Effects  of  Differences  in  Thermal  Properties  . 

4.3.3  Shearing  Stress  at  The  Coating/Substrate  Interface 


List  of  Figures 


Figure 

1.  Two-dimensional  asperity 


2. 


Thermal  stress  (ct.*)  vs.  mechanical  stress 
different  traversing  speed  . . 


for 


3.  Branch  cut  selection 


4.  Strip  of  analyticity 


5.  depths  of  the  maximum  principal  thermal  stresses  as  a 
function  of  different  Peclet  numbers  . 


6.  Thermal  stress  a  and  temperature  gradient  vs  depth 
coordinate  for  aluminum  . 


7 .  Thermal  stress  <r  .  and  temperature  gradient  vs  depth 
coordinate  for  stellite  III  . 


8.  Principal  thermal  stress  vs  depth  coordinate  for  different 
Peclet  number  . 


9.  r?m  vs  laye'-  thickness  for  varying  mechanical  properties 

(substrate  properties  are  constants)  . 

10.  Maximum  principal  thermal  stresses  in  coating  layer  vs 

coating  thickness  for  different  mechanical  mismatches 
(surface  material  properties  are  constants) . 

11.  Maximum  principal  thermal  stresses  in  substrate  vs 

coating  thickness  for  different  mechanical  mismatches 
(surface  material  properties  are  constants) . 

12 .  rj  vs  coating  thickness  for  different  mismatches  in  thermal 

conductivity  (substrate  properties  are  constants)  . 


13.  Maximum  principal  thermal  stresses  in  coating  layer  \- 
coating  thickness  for  different  mismatches  in  therms, 
conductivity  (surface  material  properties  are  constants) 


14.  Maximum  principal  thermal  stresses  in  substrate  vs  coating 
thickness  for  different  mismatches  in  thermal  conductivity 
(surface  material  properties  are  constants)  .  54 


^max  VS  coatinS  thickness  for  different  mismatches  in  thermal 
capacity  (substrate  properties  are  constants)  .  55 


^max  vs  coatinS  thickness  for  different  mismatches  in  thermal 
capacity  (surface  layer  properties  are  constants)  .  56 


Maximum  principal  thermal  stresses  vs  coating  thickness  for 
different  mismatches  of  thermal  capacity  (surface  material 
properties  are  constants)  . 


Interface  shearing  stresses  (ct  )  vs  coating  thickness  for 
different  mismatches  in  thermal^conductivity  . 


r.  -  ■  -  t.  mJ  ‘  *  *-  ■» 


Nomenclature 


Asperity  characteristic  dimension,  the  half  width  of  the 
contact  area 

Specific  heat 

The  dilatational  wave  speed  of  the  surface  layer  and  the 
substrate,  respectively 

The  shear  wave  speed  of  the  surface  layer  and  the  substrate, 
respectively 

Thickness  of  the  surface  layer 
Thermal  conductivity 
Mach  number 

Load  distribution  over  the  contact  area 

Average  pressure  over  the  contact  area 

Heat  flux  distribution  through  the  contact  area 

Average  heat  flux  through  the  contact  area 

Peclet  numbers  of  the  single  material,  surface  layer  and 
the  substrate,  respectively 

Traction  over  the  contact  area  in  the  Xj_  direction 
Temperature  field 

Displacement  in  xq  and  X£  direction,  respectively 
Dimensionless  displacement  in  X]_  and  X2  direction,  respectively 
Traverse  speed  of  asperity  (x]_  direction) 

Convective  coordinates  fixed  to  the  moving  asperity 
Material  coordinates  fixed  to  the  medium 
Coefficient  of  thermal  expansion 

The  material  region:  I  for  the  coating  surface;  II  for  the 
substrate 

Dimensionless  temperature  field  (-(T-To)k/q0a) 


7 

max 

K 

3. 

l 

£ 1 1  *  £ 1 2  >  £22 

4f 

° \ 1 '^12 >a22 
°  KVa  Zr)'a  w 

p 

6.  . 

1J 


Dimensionless  depth  of  maximum  principal  thermal 
for  single  material 

Dimensionless  depth  of  maximum  principal  thermal 
for  layer  medium 

Thermal  diffusivity 

Partial  derivative  with  respect  to  x^  coordinate 
Strain  field 

Dimensionless  coordinates  (-xj/a) 

Lame's  coefficients 

Coulomb  coefficient  of  friction 

Stress  field 

Dimensionless  stress  field 
Maximum  principal  thermal  stress 
Mass  density 
Kronecker  delta 


£  ;SS 


stress 


Poisson's  ratio 


Chapter  1 


Introduction  and  Literature  Survey 

1.1  Problem  of  Thermo -mechanical  Cracking  in  Layered  Medium 

The  present  investigation  studies  the  asperity  and  the  material 
parameters  which  influence  the  thermo -mechanical  cracking  in  a  medium 
with  a  surface  coating  layer.  The  general  failure  mechanism  is  caused 
by  the  frictional  excitation  of  a  high  speed  asperity  traversing  over  a 
coated  surface.  Maximum  tensile  stress  is  considered  the  principal 
cause  of  mode  I  cracking  in  the  hard  wear  material.  Maximum  shearing 
stress  especially  in  the  coating/substrate  interface  may  very  well  lead 
to  delamination  failure  of  the  coating  layer.  The  understanding  of 
such  a  failure  process  shall  improve  the  design  of  these  devices  by 
alleviating  the  problem  of  friction  cracking. 

When  two  bodies  in  high  speed  contact  under  heavy  loads  slide 
relative  to  each  other,  the  nominal  design  pressure  between  the  mating 
surfaces  is  based  upon  the  nominal  design  contact  area.  However,  at 
high  operating  speed,  the  actual  contact  area  can  be  reduced  by  several 
orders  of  magnitude.  As  a  result,  a  low  design  pressure  may  result  in 
a  very  high  interfacial  pressure  in  the  actual  contact  zone.  Local 
high  temperature  may  occur  due  to  excessive  frictional  heating  in  the 
vicinity  of  the  contact  regions.  Cracking  may  then  happen  caused  by 
the  combined  thermal  heating  and  mechanical  load  leading  toward  wear 
and  functional  failure  of  the  device.  This  phenomenon  is  the  so-called 


"heat  checking"  or  "thermocracking"  [1].  Both  the  reduction  .  the 
performance  life  and  the  wasted  energy  are  of  concern  to  mecnanical 
designers.  For  improvement  of  the  wear  property  of  the  surface,  recent 
effort  has  been  directed  toward  surface  modification.  Research  in 
behavior  of  coated  surface  under  asperity  excitation,  hence,  has  gained 
importance  recently. 

A  general  survey  of  the  problem  of  cracking  through  the  development 
of  a  frictional  hot  spot  was  discussed  by  Burton  [2].  A  series  of 
experiments  carried  out  by  Sibley  and  Allen  [3]  showed  photographic 
evidence  of  symmetrically  moving  hot  patches  in  the  contact  zone. 

Surface  displacements  and  temperature  field  of  a  convective  elastic 
half  space  under  an  arbitrarily  distributed  fast-moving  line  heat 
source  were  obtained,  using  integral  transform  techniques,  by  Ling  and 
Mow  [4] .  The  problem  of  thermal  stresses,  Mow  and  Cheng  [5],  and 
temperature  distribution.  Ling  and  Yang  [6],  were  examined.  Two 
dimensional  models  of  heat  checking  in  the  contact  zone  of  a  face  seal 
were  presented  by  Ju  and  Huang  [7],  Two  dimensional  models  in  the 
contact  zone  of  a  thick  layered  medium  were  presented  by  Ju  and  Chen 
[8].  In  their  two-dimensional  analysis  of  a  layered  medium,  it  was 
shown  that  the  material  parameters  affect  strongly  the  stress  state 
that  would  cause  thermo -mechanical  cracking.  The  effects  of  thin 
coating  were  shown,  but  a  parametric  study  of  a  thinner  surface  layer 


was  not  considered. 


1.2  Parametric  Effects 


The  present  work  deals  with  a  general  class  of  problems  for  thermo - 
mechanical  cracking  that  is  caused  by  a  moving  asperity  traversing  at  a 
high  speed  over  a  semi -infinite  medium  with  a  thin  surface  coated 
layer.  Such  a  thin  surface  layer  is  commonly  designed  for  surface 
modification  to  improve  the  wear  property  of  the  substrate,  such  as  a 
composite,  which  is  generally  of  poor  resistance  to  friction  wear.  The 
frictional  excitation  of  the  asperity  imparts  a  mechanical  load  as  well 
as  a  thermal  load.  Through  frictional  heating  onto  the  surface 
boundary  of  the  coated  layer,  the  high  speed  traversing  heat  source 
behaves  like  a  thermal  shock.  It  is  the  combined  thermo -mechanical 
stress  state  that  leads  to  the  initiation  of  cracks  in  the  surface 
layer,  the  substrate  or  their  interface.  The  stress  states  depend  upon 
the  asperity  characteristics,  the  coating  layer  thickness,  and  the 
material  properties  of  both  the  coating  layer  and  the  substrate.  The 
analysis  of  the  thermo -mechanical  field  in  a  layered  medium  involves 
the  matching  of  the  mechanical  and  the  thermal  properties  of  both  the 
surface  layer  and  the  substrate.  The  asperity  characteristics  also 
play  an  important  role.  The  individual  influences  can  be  protracted 
from  works  on  the  parametric  analysis  of  high  speed  asperity  excitation 
over  a  single  material  by  Ju  and  Huang  [9,10,11,12],  which  was  a  fully 
three-dimensional  model.  The  parameters,  as  inferred  from  the 
mathematical  model,  are  tabulated  as  follows; 


V 

*<Xl> 

a 


(A,p)  or  (£,*/) 

P 

J*f 
a 

K 

k 

This  paper  will  expound  their  effects  on  the  stress  states  and  on 
the  locations  of  the  maximum  stress. 

1.3  General  Theory 

At  high  speed  traversing,  the  high  temperature  and  surface  yield 
due  to  the  excitation  are  subgranular.  The  plastic  wear  and  surface 
shear  for  hard  wear  material  are  demonstrated  experimentally  to  be 
restricted  to  a  very  thin  surface  layer  (Blau  [14],  and  Ruff  and  Blau 
[15]).  If  the  critical  point  of  initiation  of  thermo-mechanical 
cracking  is  at  a  point  of  depth  by  an  order  of  magnitude  larger  than 
plastic  depth,  the  base  solid  material  subjected  to  the  asperity 


friction  is  essentially  elastic.  Furthermore,  for  the  presen"  purpose, 
the  material  is  homogeneous  and  isotropic  without  local  flavs.  The 
basic  mathematical  formulas  of  thermoelasticity  describing  the  behavior 
of  continuous  media  are  the  following: 

+  (A+ji)grad  div  y  -  (3A+2/i)a  grad  T  +  py,  (1.1) 

kV2T  -  pcT,  (1.2) 


where  u  and  T  are  the  displacement  and  temperature,  (A,/i)  are  the 
Lame's  elastic  coefficients,  p  is  the  mass  density,  a  is  the 
coefficient  of  thermal  expansion,  k  is  the  thermal  diffusivity,  and  a 
dot  over  a  variable  denotes  a  time  derivative.  The  justification  of 
the  use  of  the  uncoupled  form  of  the  Fourier  equation  is  referred  to 
Boley  and  Wiener  [16].  The  effect  of  the  dynamic  term  may  result  from 
either  a  dynamic  loading  state  or  a  thermal  shock  in  which  the  rise 
time  of  temperature  is  of  the  order  of  time  of  passage  of  the  stress 
waves  in  the  material.  Duhamel  [17]  stated  that  the  dynamic  term  can 
be  disregarded  if  the  time  rate  of  change  of  temperature  is  slow 
enough.  Parkus  [18]  showed  that  the  significant  effect  from  the 
dynamic  term  can  arise  only  when  there  is  an  instantaneous  change  in 
the  surface  temperature  or  in  the  temperature  of  the  surrounding 
medium.  In  fact,  the  dynamic  effect  is  greatly  reduced  if  the 
temperature  change  occurs  in  a  very  short,  but  finite,  interval  of 
time.  This  was  confirmed  by  Danilovskaya  [19-20],  who  studied  in 
detail  the  dynamic  effect  due  to  a  thermal  shock  on  the  surface  of  a 
half-space  and  demonstrated  that  the  maximum  dynamic  stress  is  reduced 
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to  86  percent  even  for  the  extremely  short  duration  heating  of 


seconds.  The  small  effect  of  the  dynamic  term,  in  Equation  ,1.1),  will 


be  neglected  and  wind  up  a  set  of  quasi-static  uncoupled  equations. 


The  quasi-static  equations  readily  allow  us  to  apply  the  method  of 


integral  transforms, 


Chapter  2 


Problem  Statement  and  Solution  Technique 

Under  consideration  is  the  problem  of  thermo -mechanical  cracking  in 
a  coated  medium  caused  by  a  fast  moving  asperity  whose  effect  is 
separated  into  a  moving  heat  source  and  a  moving  mechanical  load  of 
combined  pressure  and  tangential  friction  force.  The  size  of  the 
asperities  are  of  the  order  of  1  mm;  the  total  thickness  of  the  medium 
including  both  the  hard  wear  coating  layer  and  the  substrate  is  at  least 
an  order  of  magnitude  larger.  Mathematically,  the  material  is 
represented  by  a  half  space  with  the  asperity  traversing  over  the 
surface  boundary  at  a  uniform  speed  (V)  as  shown  in  Figure  1. 

2.1  Basic  Equations 

Two  sets  of  coordinates  are  considered:  x^'-xo’  re  fixed  to  the 
medium  (the  material  reference  frame),  x^-X2  are  fixed  t'  oving 

load  (  the  convective  reference  frame).  The  governing  difi  -ial 

equations  are  the  thermo-elastic  Navier's  equation  and  the  uncoupleu 
Fourier  equation,  respectively,  expressed  in  the  material  coordinates, 

(VM^)aijuj +  Vjj"?  “  Vj  +  (3V2VViT/J’  (2-1} 

K.  3. .I?  -  T0,  (2.2) 

p  n 

where  3 i  -  3/3x£,  the  indices  in  the  subscripts  i,j,k  have  the  range 


—  X 2  are  fixed  to  the  moving  load 
-x '  are  fixed  to  the  medium 


■ubetrate 


1,2;  the  summation  convention  is  used  for  all  repeated  indices 
Roman  minuscules,  $  denotes  the  layered  region:  I  for  the  cc  icing 
surface,  II  for  the  substrate.  The  stress  field  Co^j}  is  related  to 
the  displacement  field  Cuj_}  through  the  thermo- elastic  Hookian  law; 


“if  -  Vrfij +  V8^  +  M5  •  (W/‘«  •  <2-3: 

where  Sjj  is  the  Kronecker  delta.  The  field  variables,  £u.j  3  and  T, 
as  found  in  Equations  (2.1)  and  (2.2),  depend  on  x^  and  t.  The 
configuration  shown  is  Figure  1  changes  as  time  varies.  Therefore, 
time  is  an  explicit  variable,  and  the  problem  is  transient.  However, 
by  taking  into  consideration  the  uniform  properties  of  the  materials, 
asperity  motions  and  loadings,  if  the  coordinates  are  fixed  to  the 
asperity  the  configuration  remains  invariant  to  time.  We  thus  justify 
the  invariant  states  of  field  variables.  In  other  word,  with  respect 
to  a  convective  reference  frame,  the  analytical  formulation  becomes 
"steady-state";  that  is,  there  is  no  explicit  time  variable.  The 
analytical  complexity  may  be  alleviated  with  the  use  of  the  convective 
coordinates  (xj^-X2).  Equations  (2.1)  and  (2.2)  become 


aijU|  +  (l-2^)3.ju|  -  S»^M*(l-2^)alluJ  +  2(l+iy<y.T0,  (2.4 


yjj  t£  -  V<3  jT^ , 


(2.5 


1/2 

where  u  is  the  Poisson's  ratio,  M  -  [V7  p  /p.  ]  is  the  Mach  number 
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of  shear  in  region  II,  and  3^  =  p^p^/p^p n;  i.e.  $1  =  p^ .ipi^ipiV 
_  =  1.  Time  is  no  longer  an  explicit  variable.  The  boundary 


and 


conditions  are: 


on  the  surface ,  X2  -  0 


e2L  -  -*i 


( 2.6 ; 


k^T1  -  -VS2 


(2.7! 


where  kj  is  the  thermal  conductivity  in  the  surface  layer,  = 

and  R2  -  9c>(x1)  in  the  contact  region  and  zero  elsewhere.  The 
convective  heat  loss  at  the  free  surface,  being  of  small  order,  is 
neglected  without  loss  of  generality.  Regular  conditions  hold  at 
infinity;  that  is,  at  infinity 


aj  -  0  ,  T0  -  0 

J-J 


The  mechanical  and  the  thermal  fields  must  also  satisfy  the 
continuity  conditions  at  the  layer/substrate  interface,  X2  -  H 


u?  -u?^  ,  cr0? 
i  l  2j 


T1  -  T11  ,  k^T1 


(2.10 


It  is  noticed  that  the  governing  differential  Equation  (2.4),  and 
the  boundary  conditions  (2.6)  and  (2.7),  are  both  non -homo gene ous .  By 
separating  the  fields  into  two  parts,  such  that 


(  o  .Q  ,  u$  3  -  (  a  .P  ,  u|  }_  +  (  a  P  ,  u@  }  , 

ij  l  ij  l  T  ij  l  M 


(2.11 


where  the  subscripts  T  and  M  represent  thermal  and  mechanic  ;  1 , 
respectively.  The  first  set  shall  satisfy  the  non-homogeneous 
differential  equations  with  the  homogeneous  boundary  conditions;  while 
the  second  set  shall  satisfy  the  homogeneous  differential  equations  but 
with  the  non-homogene.ous  boundary  conditions.  Essentially  the 
solutions  are  delineated  into  one  resulting  from  the  temperature  field 
and  another  from  the  mechanical  traction.  They  are,  respectively,  the 
thermal  stress  field  and  the  mechanical  stress  field.  Nevertheless, 
for  single  material  (equations  and  solutions  see  Appendix  A) ,  the 
thermal  stress  (cr^)  for  various  speed  and  mechanical  stress  (a**  )  are 
shown,  in  Figure  2,  that  if  the  asperity  speed  is  larger  than  0.127 
ms'-*-  (5  in/sec),  the  thermal  stress  dominates  the  failure,  and  the 
mechanical  stress  becomes  less  important.  Therefore,  for  the  case  of 
"hot  spot",  because  the  asperity  traversing  speed  is  much  higher  (more 
than  5  ms'*-)  than  the  indicated  limit  speed  (0.127  ms"*-),  the  study  in 
the  sequel  will  focus  essentially  upon  the  thermal  stress. 


2 . 2  Governing  Equation 


By  using  the  following  dimensionless  variables: 


i  -  x^a,  rj  -  Xj/a,  D  -  H/a,  u  -  u^a,  v  -  u2/a, 


a jl'  ”  <7n/*>0  ’  '  arjrl  ~  a22/*’o> 


Cx  -  [(An+2MII)/pII]V2i  c2  -  On/Pn)1/2.  C*  -  [(Ai+2mi)/pi]1/2| 
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Figure  2.  Thermal  stress  vs.  Mechanical  stress  for  different  traversing 
velocities . 


C*  -  (A j/Pj)l/2.  b  -  b*  -  [(3Ai+2Mi\  .^l1/2, 

7 1  —  ^2  ”  *lo I  ’  ^  ”  ^1/^2 >  ^  “  ^2/^2  >  ^  > 

M  -  V/C2,  N  -  C^,  P  -  P/P0,  Q  -  q/q0,  Rj  -Va/Kj,  -Va/*.^, 

^  -  (T^-T0)kI/q0a. 

Equation  (2.4)  becomes: 


d2u£ 


32V^  ,  a2u^  .  (2.12-a) 


^  3£2  +  (a0-b£)  3^3r?  +  3r?2  Cp  d£  -  ~/9“  3^2 


,  a2v/^  a2u<° 

P  dt2  +  ^P'  p )  3f3»7 


-4-  3.  _ - 

P  dr] 2 


_  SJ1*  -  <2.i: 


32v^ 


^  .  _8^  3  „2 _ 

P  dn  ~  P  8£2 


VVvIj)  ' 

arWV  ■ 

(n2 , 1  ,b 

/c;,lj .  for  0  -  II, 

where  M  is  the  Mach  number  in  the  substrate,  and  equal  to  V/C2 ,  C2 

is  the  shear  wave  speed,  and  V  is  the  traversing  speed  (or  the  asperity- 

speed)  ,  which  is,  without  loss  of  generality,  of  the  order  of  15  ms'-'-. 
1/2 

M3  denotes,  therefore,  the  Mach  number  in  the  surface  layer,  which 

for  hard  wear  material,  such  as  Stellite  III,  is  of  the  order  of  10'^. 

2 

Since  the  magnitude  of  3M  is  of  small  order  to  a^  and  b^,  the  dynamic 
terms  on  the  right  side  of  equal  sign  in  (2-12-a)  and  (2-12-b)  are 


negligible  comparing  to  the  first  terms  of  both  equations 

Equation  (2-12)  thus  is  simplified  to  the  quasi-static  form: 


+ 


hP~!Tr  +  ( W~a?ar  + 

The  Hookian  law,  Equation  (2.3),  is  thus 


(2.13-b) 


a  n2  r  su/3  av/3 

“  TT  V  aP~W  +  (VV“*T  '  V 


(2.14-a) 


ft  m2  r  3v£  3ufi  1 

%  “  Po  W  3^  +  a?  J’ 


(2.14-b) 


ft  m2  r  duP  3v£  I 

“  “p^  V  (VV  3?  +  *0  dr,  '  V^J 


(2.14-c) 


where  d^-p^/p  that  is,  d^  -  1;  -  5.  The  boundary  conditions  are: 


-  MfP<0  at  r,  -  0,  -1  <  C  <  1, 


|  -P(0  -1<|<1  r?  -  0 

1  0  -1  >  £  and  f  >  1  rj  -  0 


(2.15-a) 


(2.15-b) 


aJJ,  aJP,  ct^P  -»  0  as  £2+r]2  ®, 

ZZ  Z.ri  VV 


(2.16) 


crJ-  -  pJ1,  a  1  -  a11  at  ij  -  D, 

€>?  nn  nn 


( 2 . 17  -  a) 


VP  -  u11,  -  vIP  at  r)  -  D. 


(2.17-b) 


The  Fourier  equation,  Equation  (2.5)  is  expressed  in  the 
dimensionless  form  as: 


In  the  next  section,  the  solution  of  these  sets  of  equations, 
Equations  (2.13),  (2.14),  (2.18)  will  be  obtained,  using  the  Fourier 


transform  method. 

2.3  Integral  Transform  Technique  and  General  Solution 

In  this  section,  we  will  develop  in  the  Fourier  transform  space  the 
general  solutions  for  the  quantities  ,  v£,  a  A,  a  P  o 

CC  C'?  nv 

2.3.1  Fourier  Transform 

In  view  of  the  boundary  conditions  and  the  problem  itself,  by  the 
results  in  [21],  the  appropriate  representation  of  the  solutions  may  be 
obtained  by  the  Fourier  transform  defined  by  the  following  equations: 

p+co 

g(s , *7 )  -J  g(C.»7)exp(tsOd£  (2.19-a) 

-00 

1  r+00+  4  C  - 

g(C.n) - 27r~  J  g(s , »?)exp( -  4s£)ds  (2.19-b) 

-  <*H-  4  C 

where  s  is  the  complex  transform  parameter,  and  g(s,rj)  denotes  the 
Fourier  transform  of  g (£,77).  Equation  (2.28)  forms  the  Fourier 
transform  pair. 

2.3.2  Mechanical  Stress  Field 
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Aik 


■  •** »  «  1 m*nJ  "h  ""*■  **  4  • 


sa 


i 

VJ 


I 


With  reference  to  Equations  (2.11),  (2.13),  and  (2.14),  ~  .= 
equations  to  be  used  for  the  mechanical  portion  of  the  solutions  are 


the  following: 


32u£  .  32v£  ,  32u£ 

+  +  h0~8^-  “ 


^~de~ +  (VV' 


d^dr)  +  &  dr)2 


0  42  .  r  <*u 

a55  “  p0  V  35 


r  dufi  3v^  -I 

L  aj3  30  +  (aB+hBy'  3n  \  ’ 


0  42  ,  ,  f  du/3  -I 

“  Po  W  35  +  ar?  J1 


a  r  3u^  3v^ 

%  "  "PT  V  (W~3T  +  a/3~ 3rT 


(2.20-a) 


(2.20-b) 


(2.20-c) 


(2.20-d) 


(2.20-e) 


The  boundary  conditions,  the  regularity  conditions,  and  the  continuity 
conditions  are,  correspondingly, 


o ^  -  /ifP(5)  at  rj  -  0,  -1  <  5  <  1, 


/  -P<0  -1<(<1 

W  1  0  -1  >  5  and  f  > 


V  “  0 


5  and  5  >  1  r?  -  0 


cr££,  <7^  -»  0  as  52+rj2  -»  ®, 

55  s'?  W 


C7  1  -  O11 


ul  -  UII, 


at  rj  -  D, 


at  r)  —  D. 


(2.21-a) 


(2.21-b) 


(2.22) 


(2.23-a) 


(2.23-b) 


vfr.V  vXv  »*>  ■- 


X-:.vX-X  X- 


9 


a 


The  transformed  equations  may  be  obtained  by  taking  the  -tie] 
transform  of  Equations  (2.20-a)  through  (2.20-e)  with  respe:~  to  £ : 


d2u^  ,  dv$  -a 

V*p-  -  <<*)<*/!■  V—  '  B  ^  '  °' 


d2v^  ,  du^ 

Vs?-  ■  «*><VV—  •  V  ^ 


-  -T7  d^[  (V2V~5r  ■  (is)V*/3]' 


P  pL  d. 


-  fl  M2  r  avp  .  1 

‘  —  V  V3T  +  <V2V^  1 1 


The  boundary  conditions  in  the  transform  space  are 


(2.24-a) 


(2.24-b) 


(2.24-c) 


(2.24-d) 


(2.24-e) 


<7  1  -  /ifP(s)  -  J4"”  P*(^)exp(  ts^)d$  at  »j  -  0,  (2.25-a) 


a  1  =.  -PCs)  - 
VV 


-  J*"  P*<Oexp(  ts?)d?  at  rj  -  0,  (2.25-b) 


where 


f  *P(0  -1  <  f  < 

p*(o  - 

l  0  -1  >  £  a 


1  >  £  and  £  >  1 


Because  the  Fourier  transform  requires  the  function  in  the  £  - 
direction  to  be  either  exponentially  bounded  or  to  vanish  at  infinity, 
therefore,  the  regularity  conditions  in  the  transform  space  at  infinity 
become : 


INV.V.WT 


'Cd 


Vj 
V ( 


O 


;ii  -  o 

C?  nv 


as  r)j 


(2.26) 


The  boundary  conditions  become: 


a  I  -  o11 

in  in' 


ul  -  ull, 


(7  I  -  ^11 

nn  m 


at  tj  -  D, 


,1  _  VII  at  i|  -  D. 


(2.27-a) 


(2.27-b) 


From  Equations  (2.24-a)  and  (2.24-b)  we  have 


^  - 


d3v^  (aa-2b  ) 


Vb0 (  s“  ’  h’3 


(-T1-)' 


(a^-  b^)  s3  dr/  ’ 


(2.28) 


(2.29) 


:  -  -  2s2-  +  s*vP  -  0. 

dr/  *  dr/  ^ 


The  general  solution  for  Equation  (2.29)  is 


-  a*  exp  ( -  r/s )  +  b*  exp(r/s)  +  c*r/  exp(-r/s)  +  d*r/  exp(r/s)  .  (2.30) 


Substituting  Equation  (2.30)  into  (2.28)  we  have 


uP  -  i  \  a*  exp(  -  r/s)  -  b*  exp  (r/s) 


+  -  D  — (  c*  exp(  -  r/s)  + 

P£  s  p 


+  d*  exp(rjs)  +r/ (c  *r/  exp(-r/s)  -  d*r/  exp(r/s)]j,  (2.31) 


where  p  -  1-a  /b  .  Then,  substituting  Equation  (2.30)  and  (2.31)  into 
P  P  P 

(2.24-c)  through  (2.24-e),  we  obtain  the  mechanical  stress  field  solution 
for  Re[s]>0,  as  follows: 


I 
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0^| - d^s[  2b p  (a*  exp(-rjs)  -  b*  exp(rjs)} 


+ 


2a.-b  1 

2b^(  b  .  a  CJ  exP<_,»a>  +  d*  exp(rjs)]  + 

PP 


+  2b^r;(c*  exp(  -rfs)  -  d*  exp(r;s))  ,  (2.32) 


P 


- ^T~  (‘s)d^b^[  -  2(a^  exp(  -T\s)  +  b*  exp(rjs)] 


P  P 


P 


-  2((  — - 1)-T-  +  n)c*  exp(-fjs)  + 

P  p  P 


+  -1'>~  -  r>)d*R  exp(ijs)]], 


P 


-  p 

a  —  - 


n  p0  p  p" 


d^b^s[  -2  (a*  exp(-f?s)  -  b*  exp(rjs)] 


—  (c*  exp(-r;s)  +  d*  exp(f7S>] 


fP 


-  2r\  (c*  exp(-r?s)  -  d*  exp(r;s)]  . 


Using  the  regularity  condition  (2.26) 


b*  -  d*  -  0. 


II  II 


(2.33; 


(2_34; 


(2.35 


For  Re[s]  <  0  we  can  get  a  similar  set  of  solutions,  where  a*,  b*, 
and  d*  are  func 
profile  P  given  by 


c*,  and  d*  are  function  of  s  and  depend  upon  the  pressure  distribution 


2sin(s) 


(2.36) 


for  a  uniform  unit  pressure  in  the  contact  region,  -1  <  f  <  1. 


Now  enforcing  the  boundary  conditions  (2.25)  and  (2.27),  the  six 
coefficients  a*,  b*,  c*,  d*,  a*^ ,  and  c*^  can  be  solved  from  six 
algebraic  equations.  The  six  algebraic  equations  may  be  presented  in 
matrix  form  as  Equation  (2.37),  and  to  be  solved  in  Chapter  3  by  Gaussian 
elimination  using  the  computer  subroutine  "DECOMP". 


Let  Ap  =  C(2/p^-l)/s4-D]e-Ds 

=  C(l-l/p^)/s-D]e-Ds 

*0  =  -C(l/(p^s>+D]e-Ds 


=  C(2/pI-l)/s-D]e-Ds 
®  =  C(l-l/Pl)/s+D]e-Ds 
7  =  -C(l/(p  s>-D]e-°s 


and 


A  UM  - 


(2.37) 


where 


A 


-1 

-1 

(l-l/Pl)/s 

-(l-l/Pl)/s 

0 

0 

-1 

1 

-1/PjS 

-1/PjS 

0 

0 

e-Ds 

-  eDs 

*1 

-e'Ds 

e*Ds 

eDs 

De‘Ds 

DeDs 

-e^s 

De  "Ds 

e’Ds 

.eDs 

e'Ds/SJ2 

-«TT/SJ2 

2.3.3  Temperature  Field 


The  governing  differential  equation  for  the  temperature  field, 
Equation  (2.18),  is  restated  here: 


d2<t>P  d2<j>0  d<j>P 


d<Z2  '  3ri: 


P  3? 


0-1,  II. 


(2.18) 


The  boundary  conditions,  the  regularity  conditions,  and  the 
continuity  conditions  are 


3$? 


dv 


f  -Q(?) 

Qx(0  - 

l  o 


-l  <  i  <  i 
- 1  >  £  and  f  >  1 


0 

o’ 


(2.38) 


d<j>‘ 


d<t> 


II 


dr) 


3 


dr) 


n  ~  D, 


(2.39-a) 


-  4>11, 


ri  -  D, 


(2. 39-b) 


I  ,*11 


4>l ,  <t> 


0,  as  i2+rt2 


(2.90) 


where  B  -  k  /k^ .  Equation  (2.18)  in  transform  space  becomes 
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■  The  solution  of  Equation  (2.50)  is 

1 

V  ^  “  A^exp(-7s2- iR^s  ?7)  +  B^exp(,/s2  -  <:R^s  17 ) 

(2.42) 

K  where  Re[,/s2-tR  s]  >0,  A  and  Bd  are  function  of  s  and 

^  PI 3/3 

&  upon  the  heat  input  profile  Q*  given  by 

depend 

I 

.  2sin(s) 

>  Q  - 

(2.43) 

£ 

P  uniform  unit  pressure  and  friction.  Using  the  regularity 

S'  condition  (2.40) 

* 

i  Bn  - 

S 

(2.44) 

S' 

,v 

Applying  the  boundary  conditions  (2.38)  and  (2.39), 

1  coefficients  Aj ,  Bj  and  Ajj  can  be  readily  found: 

■/ 

the  three 

1  ..  Q*(FI+SFn) 

£  1  Fj(l-exp(-2FID)j  +  SFjF^ (l+exp ( - 2FjD) )  ’ 

F» 

(2.45-a) 

r> 

£  n,  Q*(FI+SFn)exp(-2FID) 

*  1  Fj(l-exp(-2FID)]  +  SFIFII[l+exp(-2FID)}  ’ 

V* 

(2.45-b) 

2Q*exp((Fn-FI)D) 

&  .  11  Fj(l-exp(-2FID))  +  BF^ (l+exp( -2F  D))  ’ 

(2.45-c) 

i>fci 


where  -  Js2-i R^s,  F^  -  ,/s2-tRjjS,  Re[F^  ]  >  0,  and  Re’  >  0. 

2.3.4.  Thermal  Stress  Field 


From  Equations  (2.11),  (2.13),  and  (2.14),  the  equations  for  the 


thermal  component  of  the  solutions  are: 


32u#  ,  32v£  ,  32u£  84>P 

a3^7T-  +  +  ?~d^~  "  C^“ 


,  32v/3  32u^  d2vP  3<i 

rW~  +  (a/3'b/3)~a^“  +  ai3~J^r  -  c0—t 


3  ^  J  r  <*u 

“  P0  -3 L  a3  3? 


dufi  3v/3  .-i 

sF"  +  (£VV  3 17  '  J  ’ 


3  ^2  f  3v^  5u/3  -I 

in  P0  d/3b0i-  +  °n  J  ’ 


3  42  r  3u^  3v^  1 

°  nn  '  P0  d/3  L  (VV  3?  +  a/3~FT  *  V  J 


(2.46-a) 


(2.46-b) 


(2.46-c) 


(2.46-d) 


(2.46-e) 


The  boundary  conditions,  the  regularity  conditions,  and  the  continuity 


conditions  are 


<7.^-0  atrj-O, 
t  n 


a^-0  atr?-0, 

nn 


■  1  <  f  <  1 , 


-1  <  f  <  1 
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(2.47-a) 


(2.47-b) 


(2.48) 


(2.49-a) 


at  r\  -  D. 


2.49-b) 
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UI-UH.  V1  -  V11 


Equation  (2.46)  in  the  transformed  space  become: 


,  d2u£  ,  dv/3 


av^  o  a  - 0 

_  -  a.pS~uP  _  .(is)c^, 


(2.50-a) 


U-  vr  uur  9  « 

l/Td?“  '  -  V  "  -  <7 


(2.50-b) 


~  3  4 2  f  dvd  -o  -  I 

■’«  ‘  ~  V  (V2V^T  '  <is)V^  •  Vr 


(2.50-c) 


-  3  "2  ,  .  r 

%  "  P0  W  dr; 


-  ( <s)v/3j  , 


(2.50-d) 


-  B  r  dvP  ~a 

%n  '  —  V  VST  *  <V2V^  ■  V 


(2.50-e) 


The  boundary  conditions,  the  regularity  conditions,  and  the 
continuity  conditions  are  transformed  into  the  expressions  (see 
comments  before  equation  (2.26)) 


(2.51-a) 


a ^  -*  0  as  r}‘ 

ZZ  Zn  w 


(2.51-b) 
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(2.53-a) 


7  -  ull. 


7  —  '/II  at  rj  —  D. 


(2.53-b) 
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From  Equations  (2.50-a)  and  (2.50-b)  we  have 


if  -  - 


Vb* 


i  d3v£  (aa-2b0)  <  dv/3 

^  s3  '  d»7 3  (a^-  (  s3  '  dr? 


b*c*  1  a2 ^  ca  1 
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(2.54) 


d^v^  d2v^  _  d3<j>P  d<j>P 

-  2s2  -  ,  +  s*vP  -  E„ - E„s2- 


dr;4  dr; 2 


~P  dr)z  “/T  dr]  ’ 


(2.55) 


where  E„  -  ca/a.a.  Equation  (2.55)  is  a  fourth  order  nonhomogeneous  ordinary 
P  P  P 

differential  equation,  where  the  nonhomogeneous  part  comes  from  temperature 
field.  The  complementary  cwP  solution  of  Equation  (2.55)  is 


cv^  -  A*  exp (-tjs)  +  B*  exp(rjs)  +  C*rj  exp(-rjs)  +  D*y  exp(r;s)  .  (2.56) 

P  P  P  P 


Let  the  particular  solution  Pv/3  of  Equation  (2.55)  be  given  by 


Pv/3  -  'Tip  exp(-F^s)  +  0 p  exp(FpS)  . 


(2.57) 


Substituting  Equation  (2.57)  into  (2.55),  we  obtain: 


■2e 


•f  f;-s=  "p 


Aa  ,  0 


Yp  b' 

F=-sF  V 


(2.58) 


The  general  solution  of  Equation  (2.55)  is  obtained  by  combining  the 
complementary  and  the  particular  solutions.  That  is: 
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v.v. . 

■  vaA 


VP  -  A*  exp (-rjs)  +  B*  exp(Tjs)  +  C*r?  exp(-fjs)  +  D*r?  expCrjs 


+  *  exp(-F  s)  +  n  exp(F.s). 


Substituting  Equation  (2.59)  into  (2.54)  we  have 


—  r  2  1 

U|5  "  ‘L  Afl  exP(*»»s)  -  exp(r?s)  +  ( -  -  1) - (  C* 

r  P  P  a  s  P 


(2.59) 


exp(-»7s)  + 


+  D*  exp(rjs)  +r)[c*n  exp (-r)s)  -  exp(»js))  - 
-  Gp  exp(-F^s)  -  Hg  exp(F  s)j , 


(2.60) 


where 


-  [rf  + •, ♦  <-^>2  •  i]  ■ 

H0  -  hr  +  l)hN  % +  *,[4-  h^2  ■  l]  4-  • 


Equations  (2.59)  and  (2.60)  together  with(51-c) , (51-d) ,  and 


(51-e)  yield 


CT4  “  P*  Vl  2Va/3  exP(^s>  *  b0  exP^s>]  + 


![  2bJ; 


2a  -b  1 

+  2b/3  (  b  -  a  c3  exP(-’?s)  +  d*  exp(»7s))  + 

P  P  P  p 


+  2hpn^C*p  ^P*'*®)  '  exp(r?s))  + 


+  exp(-F^s)  -  Lg  exp(F^s)] , 


(2.61) 
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Using  the  regularity  condition  (2.61) 

bII  “  dII  “  °‘  (2,64) 

Applying  the  boundary  conditions  (2.51)  and  continuity  conditions 
(2.53),  the  six  coefficients  A*,  B^,  C* ,  D*,  A^  ^ ,  and  C*^  can  be  solved 
from  six  algebraic  equations.  They  have  been  represented  in  matrix 
form  as  equation  (2.65),  to  be  solved  in  Chapter  3  by  Gaussian 
elimination  using  the  computer  subroutine  "DECOMP". 

A  UT  -  5X  (2.65) 

where 

Transport  matrix  of  -  {  A*  B*  C*  D*  AT*  CT*  } 

and 


Chapter  3 


Numerical  Solution 

3 . 1  Selection  of  Branch  Cuts 

In  discussing  a  correspondence  v  =  A(z>  we  use  the  word  "function" 
to  mean  that  A  assigns  a  single  value  u/  to  each  permissible  value  of  z. 
Sometimes  this  fact  is  emphasized  by  saying  "  A  is  a  single-valued 
function".  Of  course  there  are  equations  which  do  not  define  single¬ 
valued  functions.  For  example,  in  the  solution  of  temperature  field 
(2.51),  Fp  =  Js2-iRps  is  not  a  single-valued  function  in  the  complex 
s -plane.  Indeed,  for  each  nonzero  s,  there  are  two  distinct  values  of 
Fg.  Since  the  concept  of  analyticity  was  defined  only  for  single -valued 
functions,  we  can  not  directly  discuss  this  property  for  multiple-valued 
fucntions  such  as  Fp  =  Js2-iRps.  We  shall,  however,  investigate  the 
analyticity  of  certain  single-valued  functions  which  can  be  derived  from 
a  multiple-valued  function  w  —  A(z).  We  construct  these  single-valued 
tunctions  by  focusing  our  attention  on  a  domain  D  and  selecting  a  single 
value  function  which  is  analytic  in  D,  and  then  we  call  it  a  BRANCH  of 
h( z)  in  D. 

To  define  the  branch  of  the  double -valued  function  on  Fp  =  J s2-  i-Rps 
such  that  Re[F^]  >  0 ,  we  write 

^  -  (s2-  ;R^s)1/2  -  S1/2  (s- tR^s)1/2,  (3.1) 

and  carefully  define,  using  Figure  3,  s1/2  -  l1  exp (t^),  (s-t'R^)1/2  - 
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l2  exp (i&2)>  wich  the  values  for  and  d2  shown  below.  On  an 
select  the  branch  cuts  along  the  imaginary  axis  in  the  s--  Lane,  and 
the  range  of  9lt  92  such  that 

-n/2  <  6X  <  2n/2 

and 

-  3ir/2  <  02  <  n/2 
so 

F0  -  (t1-£2;1/2exp[4(51+52)/2) 

Note:  1.  For  the  points  along  the  real  axis  and  along  the 

imaginary  axis  between  (0,R^),  -rr/4  <  arg(F^)  <  0. 

2.  For  the  points  just  left  to  the  cut  (0,-®),  or  left  to 
the  cut  (R/j,00) ,  arg(F^)  -  n/2 . 

3.  For  the  points  just  left  to  the  cut  (0,-®),  or  left  to 
the  cut  (R^,®),  arg(F^)  -  -n/2. 

3 . 2  Outline  of  the  Computer  Program 


(3.2-a) 

(3.2-b) 

(3.3) 


In  simple  applications  of  Fourier  transforms,  the  transform 
parameter  s  is  a  real  variable  and  the  path  of  inverse  Fourier 
transform  coincides  with  the  real  axis.  In  view  of  the  expression  in 
Equations  (2.42),  (2.61),  (2.62),  and  (2.63),  there  are  branch  point  at 


the  origin  an  at  s  -  tR^.  There  is  also  a  pole  at  the  origin.  The 


integral  is  therefore  singular.  According  to  the  branch  cut  selection, 
we  can  have  a  strip  of  analyticity  as  shown  in  Figure  4.  Knowing  from 
complex  analysis  [21]  that,  if  the  function  g(s, ri)  is  analytic  in  a 


Figure  4.  Strip  of  Analyticity 


domain  D,  the  integral  is  independent  of  the  path.  Therefor-  ichin 
this  strip  of  analyticity,  we  can  translate  the  integral  p:\  parallel 
to  the  real  axis  over  a  distance  c  to  avoid  the  singularities. 

All  solutions,  obtained  in  Chapter  2,  are  in  the  transform  space. 
Physical  solutions  are  obtained  with  inverse  transforms.  In  general, 
che  general  solutions,  Equations  (2.30,31,32,33,34,42,59,60,61,62,63), 
are  too  difficult  to  invert  analytically.  Numerical  inversion  is 
possible  using  the  computer  subroutines  "DECOMP"  to  perform  Gaussian 
elimination  and  "KQUAD"  to  perform  the  inverse  integration.  Both 
subroutines  "DECOMP"  and  "KQUAD"  are  in  a  IBM-3081D  computer.  DECOMP 
solves  a  linear  system  of  equations  by  Gaussian  elimination;  KQUAD 
integrates  real  functions  of  one  variable  over  a  finite  interval,  using 
an  adaptive  quadrature  scheme  based  upon  Gauss -Kronrod  algorithms. 

These  two  subroutines  all  deal  with  real  functions  and  variables, 
but  the  present  problem  involves  complex  functions  of  the  complex 
variable  s.  Therefore,  we  do  the  following  work  for  the  computer: 

(i)  For  the  program  "DECOMP": 

Basically,  we  solve  a  linear  system  of  equations  in  matrix  form  as 

a  s  -  e  (3.4) 


where  a,  £,  and  S  are  all  complex,  we  can  decompose  them  as 


follows 


(a  +  <b)(x  +  ty)-(c  +  (d) 


(ax  ~  by)  +  ( (bx  +  ay)  =  (c  +  (d) 


(3.5) 


This  equation  may  be  written  as 


(  a  — b"|  /  x  \  _  /  c  x 

lb  a.  ^  y  ^  ^  d 


(3.6) 


Equation  (3.6)  suggests  that,  instead  of  inputting  a  6x6  complex 
number  matrix,  we  input  a  12x12  real  number  matrix. 

(ii)  For  the  program  "KQUAD" 

We  integrate  functions  like: 


g(€.n) 


jS-oo+ 


£C 


g( 


-oo+  £C 


s,r;)exp(-  ts^)ds, 


(3.7) 


where  g is  a  real  function,  but  g ( s , »7 >  is  a  complex  function, 
Therefore,  for  a  real  variable  solution,  eventually,  we  do  the  integration 


««.•>  - «.[  p,c  g(s  ,  rj)  exp(  -  ts£)dsj  ,  (3.8) 

-oo+  ic 

where  Re [  ]  means  "take  the  real  part  of  the  function  in  the  square 
brackets",  Equation  (3.8)  can  also  be  rewritten  in  terms  of  the  real 
variable  u  as 


g  (€.»?) 


-  r  Retg(  u-  ic ,  r?)exp(-  t(u-  £C)^]]  du. 


(3.9) 


The  Equation  (3.9)  indicates  that,  instead  of  integrating  a  complex 


34 


<Vvv, 


function,  we  can  just  take  the  real  part  of  the  integrand  and  ..ms late 
tne  integral  path  by  a  distance  c  parallel  to  the  real  axis  :o  avoid  the 
singularity  at  the  origin,  and  then  integrate  it  by  computer. 


Chpater  4 


Parametric  Study 

4.1  Asperity  Parameters  and  Critical  Depth 

The  asperity  parameters  consist  of  asperity  speed  V,  pressure  P, 
and  contact  area  a.  The  coefficient  of  Coulomb  friction  /jf  relates  the 
normal  force  and  the  friction  force.  The  friction  force  generates  both 
the  mechanical  portion  and  the  thermal  portion  of  the  stress  state. 

Its  rubbing  speed  (V)  influences  the  thermal  input.  It  is  readily 
conceivable  that,  at  low  rubbing  speed,  the  mechanical  portion  of  the 
stress  dominates.  The  static  case  of  V-0  is  indeed  the  limiting  case. 
Therefore,  at  high  speed,  the  thermal  stress  prevails  (see  Figure  2). 
Without  loss  of  generality,  the  distribution  of  asperity  pressure  is 
assumed  to  be  UNIFORM,  for  mathematical  simplicity.  For  the  effect  of 
non-uniform  pressure  distribution,  see  the  results  of  Ju  and  Huang 
[9,10],  it  was  shown  that  the  assumption  of  uniform  distribution  of 
asperity  pressure  will  result  in  stresses  40%  less  than  those  from 
three-dimensional  non-uniform  distribution  of  pressure  [10,12,13], 
However,  it  is  adequate  to  use  this  assumption  in  the  study  of 
interaction  among  various  parameters  without  actually  predicting  the 
stress  level  for  fracture  initiation. 

We  shall  designate  tj  as  the  critical  depth  where  maximum 
principal  thermal  tensile  stress  occurs  for  single  material  or  for  the 
coating  as  a  single  material,  and  rj  _  that  for  the  coated  medium. 


For  various  parameters  and  coating  thickness,  the  anal-  , 
principally  will  locate  the  depth  ,?max  ■  This  depth  is  dep  lent  on  the 

material  and  the  coating  thickness  .  From  the  set  Cn  j  the  worst 

max  ’ 

case  corresponding  to  one  with  the  highest  maximum  thermal  tensile 
stress  is  important  to  the  designer.  Such  a  worst  case  is  related  to 
the  critical  depth  n 

cr 

The  depth  for  a  single  wear  material,  where  maximum  principal 

thermal  tensile  stress  occurs  is  shown  in  Figure  5.  In  this  figure, 
the  solid  line  represents  the  theoretical  t?cj_  as  a  function  of  Peclet 
number.  The  analytical  computation  of  this  solid  line  is  from  the 
solution  of  thermal  stresses  for  a  single  material  (see  Appendix  A) . 

For  the  numerical  result,  we  select  two  materials,  Stellite  III  and 
Aluminum  Oxide,  as  the  test  materials.  Varying  the  Peclet  number 
(R-Va/fc )  for  both  is  therefore  accomplished  with  varying  the  asperity 
traversing  speed.  The  principal  thermal  stresses  in  the  trailing  edge 
of  the  asperity  is  calculated  for  changing  depths  from  the  wear 
surrace.  Then  the  critical  depth  ncr  is  obtained  where  the  maximum 
principal  thermal  tensile  stress  occurs.  Hence,  we  have  the 
relationship  between  the  V cr  and  the  Peclet  number  for  the  two- 
dimensional  case. 

The  critical  depths  for  three  other  materials  Aluminum  (Al) , 
Silicon  Carbon  (SiC)  ,  and  Zirconium  (Zr)  are  computed  for  the  same 
asperity  speed  of  15  ms'1,  and  the  same  asperity  width  of  0.254  mm. 

For  each  of  these  materials,  the  material  properties  are  different  from 
those  of  uhe  test  materials.  It  is  shown  in  Figure  5,  the  critical 


fhe  Maximum  Principal  Thermal 


depth  of  each  material  falls  on  the  theoretical  lin  This 

illustrates  that  the  critical  depth,  for  a  single  asp  ,  _ty  excitation 
over  a  single  wear  material,  depends  predominantly  upon  the  Peclet 
number,  R-Va//c .  From  Figure  5,  the  relation  between  the  tj  and  the 
Peclet  number  for  the  two-dimensional  case  can  be  expressed  by  a  simple 
equation  as : 


2  275 

R(»l  )  -  20.4368 

cr 


(4.1) 


Figures  6  and  7  show,  as  functions  of  depth,  the  thermal  principal 
stress  (a'S  and  the  temperature  gradient  for  Aluminium  and  Stellite  III, 
respectively.  From  a  physical  point  of  view,  if  the  asperity  velocitv 
and  contact  width  are  constant,  higher  Peclet  number  means  the 
diffusivity  is  small.  Hence,  the  speed  of  the  heat  transfer  into  the 
medium  is  slower,  and  the  depth  rj ^  is  moderately  shallow.  On  the 
other  hand,  from  a  mathematical  point  of  view,  Equation  (A. 19)  shows 
that,  for  higher  Peclet  number,  the  decay  rates  of  the  temperature 
gradient  are  greater.  Therefore,  the  depth  to  reach  the  maximum 
principal  thermal  tensile  stress  will  approach  the  free  surface 
asymptotically  as  the  Peclet  number  increases  and  vice  versa.  Figure  6 
is  for  Aluminum  and  Figure  7  is  for  Stellite  III,  because  the  Peclet 
number  for  Stellite  III  (R  -  1400)  is  larger  than  Aluminum  (R  5  34)  for 
a  given  velocity,  consequently,  the  (rj  )  is  larger  than  (77  ) 

These  two  figures  not  only  illustrate  that  the  maximum  principal  stress 
occurs  around  the  point  where  the  temperature  gradient  begins  to 
attenuate,  but  also  imply  the  results  of  Figure  5.  The  principal 
thermal  stress  at  various  depth  at  the  trailing  edge  of  the  asperity. 
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Figure  6.  Thermal  stress  and  temperature  gradient  vs.  depth  coordinate  for 
aluminium. 


for  different  values  of  Peclet  number,  the  non-dimensional  ~  in 

Figure  8  incorporates  the  material  properties  to  yield  ma‘  . rial 
independent  values.  Therefore,  in  any  given  material,  for  a  given 
Peclet  number,  we  can  obtain  the  critical  depth  r)  from  Figure  5  or 
Equation  (4.1).  The  maximum  thermal  stress  corresponding  to  a  given 
Peclet  number  can  be  obtained  for  any  material  from  Figure  8. 

4.2  Material  Parameters 

As  the  emphasis  of  the  paper  is  on  the  friction  force  and  the  rate 
of  frictional  work,  the  simple  linear  Coulomb  law  is  used,  Fg  -  pgFn, 
for  analytical  convenience.  The  rate  of  frictional  work  heats  the 
surface  layer  of  the  wear  medium  to  raise  the  surface  temperature , 
which  in  turn  increases  the  Coulomb  coefficient  /if  with  resulting 
further  increase  in  friction  force  [22].  The  present  analysis  avoids 
the  iteration  by  using  the  maximum  Coulomb’s  coefficient  of  friction  to 
anticipate  a  steady  state  high  temperature  field.  Equation  (8)  shown 
that  the  mass  density  (p)  affects  the  stress  field  through  the  Mach 
number  (M)  and  the  thermal  field  through  the  thermal  diffusivity  («  - 
k /pc) .  Since  the  Mach  number  is  generally  so  small  for  this  class  of 
problems,  the  effect  of  the  mass  density  is  principally  reflected 
through  the  thermal  capacity  (pc),  where  c  denotes  the  specific  heat. 

The  mechanical  constitutive  coefficients  and  the  coefficient  of 
thermal  expansion  affect  the  thermal  stress  field  through  the  material 
compliances  and  the  thermal  dilatation.  Stiffer  material  and  high 
thermal  expansion  are  known  to  contribute  toward  higher  thermal 


incipal  thermal  stress  vs.  depth  coordinate  for  different 
clet  number. 


scresses.  The  effects  of  the  thermal  parameters  (k,x,pc  -  -  ^ac¬ 
tivity,  diffusivity  and  capacity)  can  be  found  from  Ju  ar.  .  Huang  [12] 
and  Ju  [13].  These  references  show  that  the  thermal  diffusivity  is 
indeed  a  derivative  thermal  parameter.  It  is  the  individual  value  of 
the  thermal  capacity,  the  thermal  conductivity  or  their  combination 
that  determines  the  thermal  field,  and  the  stress  field  that  leads  to 
thermo -mechanical  cracking  of  the  wear  material. 

Hence,  for  the  mechanical  and  thermal  impedance  matching  between 
the  surface  layer  and  the  substrate,  we  will  therefore  consider  the 
differences  in  mechanical  constitutive  coefficient,  coefficient  of 
thermal  expansion,  thermal  conductivities  and  thermal  dif fusivities . 

4 . 3  Coating  Layer  Thickness  Effects 


The  estimate  of  thermo -mechanical  cracking  from  high  speed  asperity 
excitation  in  the  layered  medium  is  more  complicated  than  that  of  a 
single  wear  material.  The  latter,  as  shown  by  Ju  and  Huang  [9,10,12], 
predicts  the  fracture  initiation  in  hard  wear  material  to  be  cohesive 
failure  at  the  depth  f?max.  where  the  maximum  principal  thermal  stress 
occurs.  For  layered  media,  the  initiation  may  occur:  (1)  inside  the 
surface  layer,  (2)  in  the  substrate,  both  through  cohesive  failure,  or 
(3)  at  the  coating/substrate  interface  through  shear  delamination.  Ju 
and  Chen  [8]  have  shown  that,  for  a  thick  coating  layer,  with  thickness 
of  the  order  of  the  asperity  size,  the  critical  stress  will  be  inside 
the  coating  at  the  thermal  layer.  The  present  paper  thus  will  address 
the  stress  states  in  the  surface  laver  as  well  as  in  the  substrate  for 
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various  coating  thicknesses,  taking  into  consideration  the  -icts  of 
mechanical  and  thermal  impedance  mismatchings .  The  folio-,  ing  results 
are  based  on  Stellite  III  as  test  material  for  both  surface  layer  and 
the  substrate. 

4.3.1  Effects  of  Differences  in  Mechanical  Impedance 

The  effects  of  both  the  elastic  coefficient  and  the  coefficient 
of  thermal  expansion  are  mechanical.  In  the  solutions  for  thermal 
stresses,  they  appear  as  a  group,  aE/(l-t/),  which  defines  the 
mechanical  impedance  for  both  the  surface  layer  and  the  substrate. 

Figure  9  illustrates  the  depths  of  maximum  principal  thermal 

stresses,  as  a  function  of  the  coating  thickness  based  upon  a  given 

substrate  property.  In  this  figure,  the  substrate  material  properties 

are  constants  and  the  mechanical  properties  of  the  surface  material  are 

considered  as  variables.  An  increase  in  value  of  the  dimensionless 

parameter  IIM  -  [aE/(l-i/)]x  /  [aE/(l-  )  ]  XX  -  0.2,  0.5,  2.0  and  5.0 

indicate,  respectively,  two  softer  surface  layers,  and  two  harder 

surface  layers.  If  coating  thickness  less  than  the  depth  rj  of  the 

coating  material  (for  our  case,  it  is  assumed  -  0.16),  the  maximum 

thermal  principal  stress  occurs  in  the  substrate.  The  rj  will 

max 

slightly  change,  depending  upon  the  value  aE/(l-i/)  of  the  surface 
layer.  If  11^  for  the  surface  is  larger  than  II  for  thesubstrate ,  then 
the  surface  layer  is  harder  than  the  substrate.  Therefore  the 
displacement  gradient  in  the  surface  layer  will  become  larger  than  the 
substrate.  From  Equation  (2.3),  indicates  that  if  the  displacement 
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gradient  of  a  surface  layer  for  a  layered  medium  is  larger  n  for  a 
single  material,  then  the  temperature  term  is  not  affecte  .  by  a  change 
in  mechanical  properties,  and  therefore  the  position  where  the  maximum 
principal  thermal  stress  occurs  will  be  shallower  than  for  a  single 
material  .  Conversely,  it  will  be  deeper  than  for  single  material.  For 
coating  thickness  larger  than  »j  ,  because  we  didn't  change  the  Peclet 
number  for  surface  laver,  the  n  will  be  the  same  as  treating  the 
surface  as  a  single  material.  Consequently,  the  layered  medium  is 
similar  to  the  single.  This  implies  that,  whether  a  layered  medium  is 
like  a  single  material  or  not,  it  depends  upon  the  Peclet  number  of 
material.  When  the  Peclet  number  is  large,  the  is  small,  and  then 

only  a  very  thin  coating  can  be  treated  as  single  material. 

In  Figure  10  and  Figure  11,  the  surface  material  is  fixed.  Figure 

10  illustrates  the  maximum  principal  thermal  stress  in  the  surface 

layer,  as  a  function  of  the  coating  thickness.  The  values  of  the 

dimensionless  parameter  IIjj  -  0.5,  1.0  and  2.0  indicate,  respectively,  a 

harder  substrate,  same  material  and  a  softer  substrate.  For  coating 

thickness  less  than  the  depth  r;  ,  the  maximum  thermal  stress  in  the 

coating  layer  occurs  in  the  neighborhood  of  the  coating/substrate 

interface.  With  increase  of  the  coating  thickness,  greater  than  the 

depth  rjcr ,  the  maximum  thermal  stress  is  in  the  neighborhood  of  the 

depth  rjmax  with  diminishing  effect  of  the  substrate.  It  is  noticed 

that,  after  a  value  of  2 r?  the  result  is  no  different  from  that  of  a 

cr  ’ 

single  material.  Figure  9  shows  the  corresponding  maximum  principal 

thermal  stresses  in  the  substrate,  a  £  ,  as  a  function  of  the  coating 

thickness.  Since  n  for  thin  coating  is  inside  the  substrate,  the 

max 
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gure  JO.  Maximum  principal  stress  in  coating  layer  vs.  coating  tliickness 
different  meclianical  mismatch. 


>vV>V 


Figure  Jl.  Maximum  principal  stress  in  substrate  for  different  mechanical 
mismatches . 


maximum  thermal  stress  will  occur  therein.  The  maximum  sub  ate 
stress  for  coating  thicker  than  the  depth  r j  will  be  at 

coating/substrate  interface.  The  prediction  of  fracture  initiation 
must  depend  on  the  cohesive  strength  of  individual  materials  of  the 
coating  and  the  substrate. 

4.3.2  Effects  of  Differences  in  Thermal  Properties 

Figure  12  shows  the  depths  of  maximum  principal  thermal  stresses, 

t]  ,  as  a  function  of  the  coating  thickness  based  upon  a  given 
max 

substrate  material  property,  with  conductivity  ratio,  -  kj/kjj,  as 

the  parameter.  The  values  of  11^  -  0.2,  0.5,  2.0  and  5.0  indicate, 

respectively,  two  more  insulating  surface  layers  up  to  two  more 

conductive  surface  layers.  In  this  figure,  we  do  not  vary  the 

substrate  material  properties.  When  the  coating  thickness  is  less  than 

rj  ,  the  value  of  tj  will  slightly  change,  due  to  the  conditions  of 
c  it  max 

the  surface  layer.  For  the  more  conductive  surface  layer,  because  more 

heat  can  be  dumped  into  the  substrate,  the  n  will  be  deeper. 
r  max 

Converserly,  r?  will  be  shallower  when  we  have  a  more  insulating 

surface.  Because  the  Peclet  number  is  not  changed,  when  the  coating 

thickness  is  larger  than  rj  ,  the  rj  will  be  equal  to  rj 

°  cr  max  cr 

Figure  13  shows  the  maximum  principal  thermal  stresses  in  the 
surface  layer,  a as  function  of  the  coating  thickness,  with  the 
conductivity  ratio,  -  kj/kjj,  as  the  parameter.  Similar  to  the 
case  of  mechanical  mismatch,  the  maximum  stresses  in  thick  coatings 
are  at  the  depth  r?  and  approach  that  of  the  single  wear  material. 


(Properties),  =  constant 


ickness  for 


When  the  substrate  is  less  conductive,  -  10.0,  the  therm  stress 
is  higher  in  the  surface  layer,  especially  on  the  coating  substrate 
-  0.5;  more  heat  is  readily  transferred  to  the  substrata.  The  thermal 
stress  is  correspondingly  reduced.  The  maximum  principal  thermal 
stresses  in  the  substrate  has  a  similar  trend,  as  shown  in  Figure  14. 

Figure  15  (substrate  material  properties  are  constant)  illustrates 

that  when  we  change  the  surface  material's  Peclet  member,  the  n  will 

be  changed  significantly.  Nevertheless,  there  is  a  general  rule  to 

help  us  to  determine  where  n  is .  We  have  used  the  rule  for  the 
r  max 

previous  two  cases.  That  is,  when  the  coating  thickness  is  less  than 

the  n  for  the  surface  layer  as  a  single  material,  the  rj  for 

layered  media  will  be  determined  by  the  substrate  material  properties, 

but  we  shall  consider  the  effects  of  surface  layer.  At  this  moment, 

when  n^c  <  1,  it  means  the  surface  temperature  gradient  is  larger  than 

the  substrate.  Therefore,  from  Equation  (2.3),  the  T7max.  i-s  slightly 

larger  than  rj which  is  defined  by  the  coating  material.  Converserly, 

the  n  _  is  smaller  than  r\  .  When  the  coating  thickness  is  larger  than 
m^x  c  it 

n  for  the  surface  as  a  single  material,  the  r  will  be  determined 
cr  max 

by  the  surface  material  properties.  Similarly,  Figure  16  shows  that, 

when  the  surface  material  properties  are  constant,  then  rj  has  the 

r  r  max 

same  tendency. 

Figure  17  illustrates  a  combined  cur/e  for  maximum  principal 
thermal  stresses,  indicating  that  the  maximum  stress  for  thin  coating 
is  inside  the  substrate,  but  it  is  in  the  surface  layer  for  thick 

In  both  cases,  the  stresses  are  evaluated  at  the  depth  rj  . 
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coat  mgs . 


Figure  14.  Maximum  principal  stress  in  substrate  vs.  coating  thickness  for 
different  mismatches  in  thermal  conductivity. 
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4.3.3  Shearing  Stress  at  the  Coating/Substrate  Interface 

Because  of  the  asymmetrical  distribution  of  the  temperature 
gradient,  the  direction  of  the  maximum  temperature  gradient  is  as 
anticipated,  not  perpendicular  to  the  wear  surface.  The  resulting 
value  of  shearing  stress  at  the  coating/substrate  interface  is  not  to 
be  ignored.  In  Figure  18,  the  curve  for  11^  -  1.0  can  be  considered  as 
the  shearing  stress  at  different  depths  from  the  wear  surface  for  a 
single  wear  material.  The  worst  case  shows  the  maximum  shearing  stress 
at  almost  30%  of  the  maximum  principal  stress.  The  value  is 
significant  in  designing  the  coating/substrate  bonding  strength  to 
avoid  shear  delamination.  The  difference  in  thermal  conductivities  of 
the  coating  and  the  substrate  is  shown  to  affect  the  interface  shearing 
stress  for  thin  coating  layers,  If  the  coating/ substrate  bonding  of  a 
thin  coating  is  a  primary  consideration,  the  substrate  may  purposely, 
by  design,  be  less  conductive. 
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Chapter  5 


Conclusions 

The  paper  demostrates  the  significant  effects  of  the  Peclet 
number  and  the  material  mismatch  between  the  surface  layer  and  the 
substrate  in  high  speed  asperity  excitation.  Principally  it  is  the 
critical  depth  rj ^  at  which  the  thermal  tensile  stress  reaches  a 
maximum.  Peclet  number  is  proven  to  dominate  the  determination  of  the 
critical  depth,  which  is  obtained  from  maximizing  the  thermal  tensile 
stress  in  the  vicinity  of  the  asperity  trailing  edge.  In  other  words, 
the  asperity  velocity,  the  asperity  size,  and  the  material  diffusivity 
are  the  principal  three  factors  that  dominate  the  critical  depth.  A 
simplified  depth  rj ^  as  Equation  (4.1)  can  be  approximated.  From  this 
result,  the  maximum  tensile  stress  in  its  alternate  non-dimensional 
form  can  be  obtained  for  all  wear  materials.  Since  the  numerical 
empirical  equation,  Equation  (4.1),  has  been  demonstrated  not  to 
depend  on  the  other  material  properties  nor  the  magnitude  of  asperity 
excitation,  the  coefficient  and  the  exponent  are  therefore  constants 
in  the  two-dimensional  model.  It  is,  however,  expected  that  the 
coefficient  and  the  exponent  in  Equation  (4.1)  would  be  function  of 
the  aspect  ratio  of  the  asperity  contact  area.  The  work  is  under 
investigation  with  a  three  dimensional  formulation. 

The  critical  depth  rj ^  is  important  in  the  layered  media  analysis. 
It  was  shown  by  Ju  and  Chen  [8]  that  the  layer  thickness  is  in  the 
neighborhood  of  rj  ,  the  thermo -mechanical  stress  state  is  most  damagin 


On  the  other  hand,  the  highest  stress  occurs  when  the  ting 
thickness  is  in  the  neighborhood  of  the  17  ,  regardless  c.  the 

differences  in  mechanical  or  thermal  properties.  As  the  coating 
chickness  increases  to  more  than  the  critical  depth  r? ^  for  the  coating 
as  asingle  material,  there  is  little  effect  of  the  substrate  on  the 
maximum  thermal  stress  in  the  coating  layer.  Hence,  if  a  thick  coating 
is  permissible  and  if  the  asperity  size  can  be  statistically 
ascertained,  a  properly  chosen  coating  layer  large  than  for  coating 

material  can  be  estimated  to  insure  the  improvement  on  wear  resistance 
for  materials  such  as  composites.  For  thermal  stress  analysis  of  a 
layered  media,  in  order  to  study  the  impedance  matching  of  the  surface 
layer  and  the  substrate,  it  is  adequate  to  consider  the  thermal 
conductivity  (k) ,  and  the  thermal  capacity  (pc)  ,  while  the  effects  of 
their  differences  are  negligible  for  thick  coatings;  the  effects  on 
thin  surface  coating  are  important.  Basically,  substrates  of  low 
mechanical  properties,  higher  thermal  conductivity  and  capacity  will 
result  in  lower  stresses  in  both  surface  layer  and  the  substrate.  The 
initiation  of  the  cracks  will  occur  wherever  the  stress  level  surpasses 
the  cohesive  strength  of  the  respective  material.  The  critical 
locations  in  the  surface  layer  or  in  the  substrate  depend  on  the 
coating  thickness  relative  to  the  *1  ■ 

The  shearing  stress  at  the  surface  layer/substrate  interface  is 
by  no  means  trivial,  depending  again  on  the  surface  coating  thickness. 
The  interface  shear  reaches  a  maximum  when  the  coating  thickness  is 
in  the  neighborhood  of  1  ■  This  knowledge  is  important  for  the  design 

of  bonding  of  the  surface  coating.  Any  shear  crack,  resulting  from  the 
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The  Equations  and  Solutions 
for  Single  Material 


In  this  Appendix  we  present  the  solutions  for  the  zeroth 
Approximation. 


(I)  Mechanical  Stress  Field: 
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The  boundary  conditions,  the  regularity  conditions  are: 
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The  solutions  are 


v  -  a,lrexp( -rjs)  +  c*r;exp  ( -  r;s )  , 
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P  - - for  uniform  unit  pressure 
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in  the  contact  region  -1  <  f  <  1 


(II)  Temperature  Field 


The  governing  equation  are: 
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The  boundary  conditions,  the  regularity  conditions  are 
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The  solution  is: 
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(III)  Thermal  Stress  Field 


The  governing  equations  are: 
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